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<^> , Abstract 

The Inert Doublet Model (IDM) , a two Higgs extension of the Standard Model with an unbroken 
Z 2 symmetry, is a simple and yet rich model of dark matter. We present a systematic analysis of the 
dark matter abundance and investigate the potentialities for direct and gamma indirect detection. We 
show that the model should be within the range of future experiments, like GLAST and ZEPLIN. The 
lightest stable scalar in the IDM is a perfect example, or archetype of a weakly interacting massive 
' particle. 
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Contemporary cosmological observations concur to indicate that the majority of matter in the universe 
not only does not shine but is not even made of ordinary atoms [H [2] . Deciphering the nature of this 
so-called Dark Matter has become one of the most important issue at the frontier of particle physics, 
astrophysics and cosmology. A profusion of dark particles have been proposed over the years and it 
is much hoped that present and forthcoming experiments will throw some light on the matter. For a 
review, see for instance OS]. 

In the present article we investigate further a rather mundane, although in our opinion quite inter- 
esting, form of dark matter. The particle candidate is a weakly interacting massive scalar and it has 
been advocated recently in |5] and [fil[7]. The framework is that of a two Higgs doublets, Hx and H 2 , 
version of the Standard Model with a Z 2 symmetry such that 
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Hx -> Hx and H 2 -» -H 2 . 

All the fields of the Standard Model are even under Z 2 . This model has been first discussed by 
Deshpande and Ma in QjJ. Following [6], we will assume that Z 2 is not spontaneously broken, i.e. 
H 2 does not develop an expectation value. Among other things, the discrete symmetry prevents the 
appearance of flavour changing neutral currents. The model is definitely not the (H 2 ) — ► limit of a 
generic two Higgs model like, for instance, the Higgs sector of the MSSM. 

The most general, albeit renormalizable, potential of the model can be written as 

V = iA\Hx\ 2 + \x\\H 2 f + \x\Hx\ 4 + A 2 |ff 2 | 4 + \ 3 \Hx\ 2 \H 2 \ 2 + M\h\h 2 \ 2 + ^ [{h\h 2 ) 2 + h.c] . (1) 

There is an Peccei-Quinn U(l) global symmetry if A5 = 0. This limit is however not favored by dark 
matter direct detection experiments (section [23)1 . 

The SU{2) x 17(1) symmetry is broken by the vacuum expectation value of Hx, 

<*> - 35 

with v = —fi 2 /Xx — 248 GeV while, assuming fi 2 > 0, 

(H 2 ) = 0. 
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The mass of the Brout-Englert-Higgs particle (h or Higgs for short) is 




2^1 = 2X 1 v 2 



(2) 



while the mass of the charged, H + , and two neutral, Hq and Aq, components of the field H 2 are given 



For appropriate quartic couplings, H or A is the lightest component of the H 2 doublet. In the absence 
of any other lighter Z 2 -odd field, either one is a candidate for dark matter. For definiteness we choose 
H . All our conclusions are unchanged if the dark matter candidate is A instead. Following [6] we 
parameterize the contribution from symmetry breaking to the mass of Hq by Xl — (A3 + A4 + As)/2, 
which is also the coupling constant between the Higgs field h and our dark matter candidate Ho. Of 
the seven parameters of the potential HJ, one is known, v, and four can be related to the mass of the 
scalar particles. In the sequel, we take [i 2 and A2 as the last two independent parameters. The later 
actually plays little role for the question of dark matter. 

The abundance of Hq as the dark matter of the universe as well as the constraints from direct 
detection experiments have been touched upon in the recent literature [HI [7] (see also [9]). In the 
present article we investigate the potentialities for indirect detection of Hq using gamma ray telescopes 
and detectors. For the sake of completeness, we also present a systematic (albeit tree-level) analysis of 
the abundance of the Hq in the light of WMAP data, as well as the prospect for direct detection by 
underground detectors. 

Before doing so, we would like to briefly emphasize some of the virtues of this, so-called, Inert 
Doublet Model (IDM), both in general and for the issue of dark matter in particular. True, the model 
is ad hoc. Also, it does not address any deep issue, say the hierarchy problem. Yet it is a very simple 
extension of the Standard Model and its phenomenology is nevertheless very rich [3 [U QUI [11]. In 
reference [6|, Barbieri et al considered this model as a way of pushing the mass of the Higgs particle 
toward the TeV scale, without contradicting LEP precision measurements (see also [10], and section [3?T1 
of the present paper). Also, in the simplest version discussed here there are no Yukawa couplings to the 
H 2 doublet, but it is somewhat natural to imagine coupling it to (odd) right-handed neutrinos. This 
opens the possibility of generating the mass of the SM neutrinos through loop corrections, a mechanism 
introduced by Ma in [5] and further addressed in a series of papers |12 [ IT3 j. IT4] . 

As a dark matter candidate, Hq is not without interest either. It belongs to the family of Weakly 
Interacting Massive Particles (WIMP) and has been proposed as a dark matter candidate by Ma . The 
most acclaimed member of this family is certainly the neutralino, the lightest supersymmetric particle 
(LSP) . The supersymmetric extensions of the SM are without doubt very interesting and well motivated 
but they involve many new parameters and their phenomenology is, to say the least, complex. It is 
perhaps far fetched to compare the respective advantages and disadvantages of the IDM and of the 
MSSM but it is interesting that the Hq shares many similarities with the LSP (weak interactions, 
similar mass scale, etc) while being considerably much simpler to analyze. Interestingly, if we take 
seriously the idea that Ho is the dominant form of dark matter, the parameter space of the IDM model 
is quite constrained. It is also, to some extent, complementary to that of the MSSM, as we show in 
the conclusions. And, cherry on the top, the phenomenology of the Ho as a candidate for dark matter 
is neatly intertwined with that of the Higgs particle. In our opinion, it is this conjunction of simplicity 
and richness that makes the lightest stable scalar a perfect example, or archetype, of dark matter. 

The plan of the paper is as follow. In the next section |(2]) we give a summary of our methodology 
for the calculations of the dark matter abundance, the photon flux for indirect detection by gamma ray 
telescopes and the cross-sections for direct detection. We then discuss in some details the results of our 
analysis, in the light of WMAP data, and study the prospects for both indirect and direct detection 
by existing and forthcoming experiments. We give our conclusions in the last section, together with 
comparison of the Model with that of the MSSM, and the prospects for future analysis. 
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2 Dark matter aspects 



2.1 Dark Matter abundance 

The abundance of Hq has been estimated in [6] and in [7] for two different mass ranges. As in these 
papers, we consider only the standard freeze-out mechanism. There are then essentially two regimes, 
depending on the mass of Hq with respect to that of the W and Z bosons. If Mh > 80GeV, the 
H annihilate essentially into Z and W boson pairs (see Figure HJ). For Mh > M h , the annihilation 
channel into h pairs opens (see Figure (J2j) ) . Otherwise, the H annihilate essentially through an 
intermediate Higgs, provided the Higgs itself is not too heavy. Some amount of coannihilation of H 
with the next-to-lightest scalar particle Aq (resp. H + ) through a Z (resp. W + ) boson can be present 
|15j . a feature that substantially complicates the determination of the dark matter abundance (see 
Figure (0)). 

We have computed the relic abundance of Ho using micrOMEGAs2.0, a new and versatile package 
for the numerical calculation of Dark Matter abundance from thermal freeze-out [16] , The latest 
implementation of this code, which was originally developed to study supersymmetric models, allows 
one to enter any model containing a discrete symmetry that guarantees the stability of the dark matter 
particle. The code takes advantage of the fact that each odd (i.e. non standard) particles of the model 
will eventually decay into the lightest odd particle, in our case the H . One can then sum the system 
of Boltzmann equations of all odd species and the relic density of dark matter is calculated by solving 
the Boltzmann equation for the lightest particle only [i"7] : 

% = \[^^M Pl (av){V\T) - r e 2 g (T)), (4) 

where Y = tidm/s is the comoving density of dark matter. All the processes of the model enter in the 
thermally averaged cross-section (ov). Integrating Eq.Q from T = oo to T = Tg the relic abundance 
is given by 

n DM h 2 = 2.72 x 10 8 ^g^r(T ). (5) 

MicrOMEGAs2.0 itself is build upon CALCHEP, a code for computing tree-level cross-sections. We 
have used both MicrOMEGAs2.0 and CALCHEP to compare with our analytical calculations and to 
perform a number of cross-checks that have strengthened our faith in the numerical results reported in 
the present paperQ 



2.2 Indirect detection 

The measurement of secondary particles coming from dark matter annihilation in the halo of the galaxy 
is a promising way of deciphering the nature of dark matter. This possibility depends however not only 
on the properties of the dark matter particle, through its annihilation cross-sections, but also on the 
astrophysical assumptions made concerning the distribution of dark matter in the halo that supposedly 
surrounds our galaxy. The galactic center (GC) region is potentially a very attractive target for indirect 
detection of dark matter, in particular through gamma rays. The produced gamma ray flux from the 
annihilation of dark matter particles can be expressed as 

*j- = y;gS <ff<t ,) \ [ P ui, ( 6) 

dQdE ^ dE^ ' inm 2 DM Ji. . s . F K ' 

where m^M is the dark matter particle mass, p is the dark matter density profile, (oiv) and dNJ y /dE 1 
are, respectively, the thermally averaged annihilation cross section times the relative velocity v and the 
differential gamma spectrum per annihilation coming from the decay of annihilation products of final 
state i. The integral is taken along the line of sight. The processes involved are shown in the Figure 
([3]), first diagram and Figures Jl]) and ([2]). 

*The Feynman rules and definition files of the IDM to be used with micrOMEGAs can be obtained upon request to the 
authors of the present article and will be included in the next available version of the code. 
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The models will be constrained by the existing EGRET [18] experimental limit on the flux, ~ 10~ 8 
photons cm _2 .s _1 , and the forthcoming GLAST [19] expected sensitivity, ~ 10~ 10 — 10~ n photons 
cm _2 .s _1 for AD, — 10~ 3 and 10" 5 srad respectively. 

Our purpose for the time being is merely to prospect the IDM with regard to gamma ray indirect 
detection. In this paper, we focus on the particle physics parameter dependence of the model and 
work within a fixed astrophysical framework, the popular Navarro-Frank- White (NFW) halo profile 
[20] . With Ro = 8.5 kpc equal to the distance from to Sun to the GC and po = 0.3 GeV.cm -3 the dark 
matter density in our neighborhood, the NFW profile density can be parameterized as 



with (a, a, (3, 7) = (20,1,3,1). It should be emphasized that the (dark) matter distribution in the 
innermost region of the galaxy is poorly known so that 7 is not very constrained, a freedom that can 
give rise to very different values of the gamma ray flux. Typically a suppression or an enhancement 
of two orders of magnitude can be obtained if one considers respectively a halo with a flat core {e.g. 
isothermal, 7 = 0) or a deeper cusp (7 ~ 1.5) |21j coming for instance from baryonic infall. Hence, 
depending on the astrophysics assumptions, the estimation of gamma ray signal can vary quite strongly 
(see e.g. [Ml HI]). 

To calculate the flux we have integrated (J6j) for a NFW profile above lGeV and around a solid 
angle of Aft = 10~ 3 srad for EGRET and Aft = 10~ 5 srad for GLAST. The differential spectra of each 
channel are given by micrOMEGAs (we have checked the agreement with Pythia simulations of |23j) 
as well as {av) at rest and we have integrated the square of the dark matter density along the line of 
sight. 

2.3 Direct detection 

A local distribution of weakly interacting dark matter could be detected [24] by measuring the energy 
deposited in a low background detector by the scattering of a dark matter particle with a nuclei of 
the detector. One distinguishes spin dependent and spin independent interactions. For Hq interacting 

2" 

with the quarks of the nuclei, there are two spin independent processes at tree level, Hoq — ► Aoq 

and Hoq — + Hoq (see Figure (@])). The experiments have reached such a level of sensitivity that the 
Z exchange contribution is excluded by the current experimental limits. Consequently, to forbid Z 
exchange by kinematics, the mass of the Aq particle must be higher than the mass of Hq by a few 
lOOkeV (we will thus disregard the A5 — > limit). We assume that the main process here is the one 
with Higgs particle exchange h. The cross section at tree-level is [6j 



where / is a form factor estimated in the literature to be / ~ 0.3 and m r is the reduced mass of the 
system. [25j. In addition, for Mh < M\y, we have neglected the one- loop exchange of two gauge 
bosons. For Mh > Mw, we include this process with, following [7j, ph -p — 4.610~ 13 pb. We take 
into account the constraints given by the CDMS [26] experiment (~ 10~ 6 pb) and the model is also 
compared to the next generation experiments like EDELWEISS II |27j or a ton-size experiment like 
ZEPLIN [28j with the valley of their sensitivities respectively around a ~ 10~ 8 and a ~ 10~ 10 pb. 

3 Analysis 

The processes driving the (co) annihilation cross section {av) relevant for the relic density from freeze- 
out and indirect detection of gamma rays are shown in the Figures {J), ([2]) and ([3]). The ifo-proton 
interactions relevant for direct detection are shown in Figure |(4]). 

There are two qualitatively distinct regimes, depending on whether the Hq is lighter than the W 
and Z and/or the Higgs boson, in which case annihilation proceeds through the diagrams of Figure 




(7) 




(8) 
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([3]). At low mass, the second diagram of ([3j) may contribute to the dark matter abundance if there is 
a substantial number of A at the time of freeze-out (i.e. coannihilation). 

In the plots of Figure ((5]), the relic abundance of Hq dark matter particles, the gamma ray flux 
due to Hq annihilation at the galactic center (indirect detection) and, finally, their cross-section for 
elastic scattering off a proton through Higgs boson (THo-p (direct detection) exchange are shown, for 
two particular Higgs masses (Mh = 120 GeV and 200 GeV) in the (Mh ,^2) plane. We refer to this as 
the low Ho mass regime. For higher Ho masses, similar plots are displayed in Figure iJH) (for Mh — 120 
GeV), respectively for the dark matter abundance, the gamma ray flux and the cross-section for direct 
detection. In each of these three cases, the color gradients correspond respectively to gradients in 
log(Q Ho h 2 ), log($ 7 (cm~ 2 s -1 )) and log(er ffo _ p (pb)). 

Since we plot our results in the (Mh i^2) plane, the diagonal line corresponds to Xl = 0, i.e. to no 
coupling between Ho and the Higgs boson. Away from this line, Xl increases, with Xl < (resp. Xl > 0) 
above (resp. below) the diagonal. Also, we write AMAq = Ma — Mh and AMH C = M H + — Mh . 

In the plots of the dark matter abundance, the areas between the two dark lines correspond to 
regions of the parameter space such that 0.094 < Q.uMh 2 < 0.129, the range of dark matter energy 
densities consistent with WMAP data. In the cross section and gamma ray flux plots, the lines indicate 
the areas of the parameter space within reach of the various experiments we take in consideration. 

3.1 Constraints 

The shaded areas in the plots of Figures J5]) and J6]) correspond to regions that are excluded by the 
following constraints: 

• Vacuum stability: 

Vacuum stability (at tree level) demands that 

Ai, 2 > 0, 

A 3 ,A 3 + A 4 - | A 5 1 > -2v/A^. (9) 

As a result, negative couplings, and Xl < among others, are largely excluded. This is the 
shaded area in the domains [ii > Mh of Figures |(5]) and J6]) . 

• Perturbativity: 

Strong couplings |Ai| > 47r are excluded but intermediate couplings, 1 < |A^| < 4ir, might 
be tolerated. The latter correspond to the regions with horizontal lines. Notice that the mass 
splitting relative to the [ii scale is smaller in the high mass regime (see Eq.j3])). For this reason, 
going away from the diagonal line in the plots of Figure ([6]) , we run faster into the large coupling 
regime than in those of Figure ^ . The regions excluded by the strong coupling constraint is of 
course symmetric with respect to the \i 2 — Mh axis. 

• Charged Higgs scalar: 

The mass of the charged scalar H + is constrained by LEP data to be larger than 79.3 GeV 
[30] . As we fix the mass differences in our plots, this constraint translates in the excluded region 
Mh < 30 GeV in the low mass regime abundance plots. 

• Electroweak Precision Tests (EWPT): 

New physics can affect electroweak precision measurements. The impact of the new H2 doublet 
can be described in term of the S, T electroweak precision parameters. As pointed out in [6], 
with appropriate mass splittings between its components, an H 2 could screen the contribution 
to the T parameter of a large Higgs masses, Mh ~ 500 GeV. The effects on the S parameter is 
parametrically smaller for the region of parameter space satisfying the previous constraints. In 
all the case we have studied, we have checked that the contribution to the T parameter resulting 
from the Higgs mass and the H2 components are consistent with EWPT. Unlike Barbieri et al, 
we only considered small values of the Higgs mass, Mh < 200 GeV. Given our choice of the Higgs 
mass Mh and of mass splittings between the H particles and the other components of H 2 , the 
constraints from EWPT are easily satisfied. 

We now turn to the analysis of the model, starting with the low Hq mass regime. 



5 



Hi 



W+ (Z) 



o , 



H 







W~ (Z) 

„2 



Hi 







^J¥+ (Z) 
H-(A ) 

^w- (Z) 

,2 



^0. 



(Z) 



W - (Z) 



9" 9" A L 5 

Figure 1: Annihilation channels into gauge bosons final state with corresponding couplings. 
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Figure 2: Annihilation channels into Higgs final state. 
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Figure 3: (Co)Annihilation channels into fermion anti-fermion final state. 
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Figure 4: Leading channels contributing to <jh -p (direct detection). 



3.2 Low mass regime, Mh around 100 GeV 

The results of this section are shown in Figure Two processes are relevant below the W, Z or 
h threshold: Hq annihilation through the Higgs and Hq coannihilation with Aq through Z exchange. 
Both give fermion-antifermion pairs, the former predominantly into bb. As the mass of Hq goes above 
W, Z or h threshold, Hq annihilation into WW, ZZ and hh become increasingly efficient, an effect 
which strongly suppresses the Hq relic density. 

Coannihilation into a Z may occur provided AMAq is not too important, roughly AM must be of 
order of Tf a ~ Mjy /25. Otherwise H annihilate dominantly into bb. This regime is quite interesting: 
the interactions are either completely known (i.e. electroweak interactions), or highly sensitive to the 
Higgs boson mass Mh and the H h effective coupling A/J 

For the sake of illustration, we study two cases, Mh — 120 GeV and Mh = 200 GeV to show how the 
predictions change as a function of the Higgs boson mass. A few general lessons can be extracted from 
these two specific examples. Note that the mass differences AMAq and AMH C are fixed respectively 
to 10 and 50 GeV while AMAq is large enough to avoid the constraints from direct detection and small 
enough to have a some amount of coannihilation into a Z . 



We may distinguish five different regions (see Figure ((Sj) , top left): 



Regions 1 and 2 are excluded respectively by charged Higgs production and vacuum stability con- 
straints. The latter involves the self coupling of the Higgs, Ai and thus depends on the Higgs mass 
(see Eq.([2])): region 2 is broader for smaller Higgs mass Mh. Although region 1 is experimentally ruled 
out, it is interesting to understand the trends. In this region, the relic abundance of Hq decreases with 
increasing Mh as one goes along [ii = Mh„ ■ Comparing the Mh = 120 GeV and Mh = 200 GeV plots 
we observe that the gradients depend, again, on the Higgs mass. This is simply due to the M^ g /M^ 
dependence of the annihilation cross section through the Higgs, i.e. smaller abundance for higher Mh 
and/or lower Mh- In region 2, the abundance is smaller the further one deviates from \i-x = Mh , re- 
flecting the dependence of the annihilation cross section on |Al|. For Mh < My/ (resp. Mn ^Mh), 
HqHq — > bb (resp. HqHq — > hh) is the dominant process. Otherwise, for Mw ^5 Mh ^5 Mh, 
H H -> WW, ZZ dominate. 



In region 3, the couplings are rather large, but nevertheless consistent with vacuum stability. Since 
\i2 > Mh q , My/, the dominant contribution to the cross section is given by HqHq — > W + W~ process, 
large enough to bring the relic abundance far below WMAP data. 



Region 4 is below the W threshold. It is the only region consistent with the dark matter abundance 
predicted by WMAP data in the low mass regime. The process that determines the relic abundance of 
Hq is again the annihilation through the Higgs. Coannihilation processes become dominant near the 
resonance at Ma +Mh ~ Mz. This is the origin of the dip in the Hq relic abundance around Mh ~ 40 
GeV (corresponding to Ma ~ 50 GeV for our choice of mass splittings) . Similarly, annihilation near 

HqH + coannihilation is suppressed for our choice of AMH C . 
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the Higgs resonance generates a second dip around M Ho w M h /2. There is an island consistent with 
WMAP data, which extends up to M Ho ~ 80 GeV, at which point WW annihilation becomes important 
and suppresses the relic abundance. 

Region 5, finally, corresponds to Mh > M w . Annihilation into a gauge boson pair is dominant, 
at least as long as Mh < M h . Gauge boson pair production dominates above the Higgs threshold for 
\i2 < Mh while Higgs pair production is dominant if ^2 > Mh , In all instances, the relic abundance 
is too suppressed to be consistent with WMAP for Mjj < 600 GeV. 

As expected, and as revealed by visual inspection, the photon flux plots shares some of the charac- 
teristics of the abundance plots. The only new salient feature is the absence of a dip at the Z resonance. 
This is of course because today, in contrast with the early universe, all Aq are gone. Existing experi- 
ments are not very constraining but future gamma ray detection experiments, such as GLAST, might 
severely challenge the model. As usual we should bear in mind the astrophysical uncertainties: by 
changing the galactic dark matter density profile, we can get higher (or smaller) fluxes. 

The plot of the cross section for direct detection is pretty transparent. Suffices to notice (see section 
I2.3[) that VHa-p at tree level is proportional to A|. The dominant contribution to the elastic scattering 
cross section is thus zero on the ^2 = Mh axis while it increases for larger values of — M/iol- 

The dependence on the Higgs mass (Eq. [8]) clearly appears when comparing the <th -p plots for 
Mh = 120 GeV and Mh = 200 GeV. For Mu > Mw, the one-loop contribution to the cross section 
from W exchange is taken into account. However it does not affect much the results as it only amounts 
for 10 -13 pb. Unless we suppose that the mass of Hq and Aq are nearly degenerate, existing direct 
detection experiments are not very constraining. Forthcoming experiments however might put a dent 
on some of the solutions consistent with WMAP. 

3.3 The high mass regime, M H 3> M w 

The results of this section are shown in Figure ©, left column. No new annihilation channel opens if 
Mh is heavier than the Higgs or gauge bosons. There are then essentially two sort of processes which 
control both the abundance and the gamma ray flux: the annihilation into two gauge bosons, dominant 
if fi2 < Mh q , and the annihilation into two Higgs, which dominates if fi2 > Mu - Coannihilation plays 
little role. It affects a bit the relic abundance along the diagonal but, even so, it is not the key process 
to get the WMAP abundance. 

The abundance of dark matter is suppressed over most of the area of the plot because of large 
quartic coupling effects on the cross sections. Strong couplings are excluded on a physical basis, but 
we found this limit nevertheless useful to understand the interplay between the various processes. In 
the present subsection, we will argue that it is possible to reach agreement with WMAP data, but only 
to the price of some fine tuning between the different annihilation channels. 

Consider first the annihilation into two Higgs bosons (the dashed line in the graphs Figure ([7])). Its 
cross section is vanishing for Al = 0. For /i 2 < Mh (corresponding to Al > 0), there is a destructive 
interference between the diagrams of Figure (J2j) , which is absent if (12 > Mh - The annihilation 
into gauge bosons depends on the quartic couplings between the scalars (see Figure ([1])). Indeed, the 
annihilation through an intermediate Higgs is controlled by A^. Also, the t and u channels exchange 
diagrams are sensitive to the mass differences between the components of the H2 doublet, which depend 
respectively on A5 for the annihilation into Z bosons and on A4 + A5 for the annihilation into W bosons. 

If Al — 0, there is no or little annihilation into a Higgs. The relevant diagrams are then the quartic 
vertex with two gauge bosons and the t and u channels with A (resp. H + ) exchange. If A5 = (resp. 
A4 + A5 = 0) the cross section into a Z (resp. W) boson pair is minimal and scales like cq^jM^ , 
If, for instance, A5 ^ the annihilation into a Z boson pair receives a contribution which grows like 
a^ (/ (M J 4 — Mh ) 2 /M%. A similar result holds for the annihilation into a W boson pair provided 
A4 + A 5 7^ 0. Notice that a^(M J 4 — M ffo ) 2 /M| ~ A|/Af^ o . The effects of weak isospin breaking 
between the components of H% is reminiscent of what happens for the Higgs in the regime of strong Ai 
coupling (the so-called Goldstone boson equivalence theorem |31j). 
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To have an abundance of dark matter in agreement with WMAP, the mass splittings between the 
components of H 2 must be kept relatively small. First because large mass splittings correspond to 
large couplings and second because the different contributions to the annihilation cross section must 
be suppressed at the same location, around \l = (ie Mh — {J< — Ma — Mh + in this case). This 
is illustrated in Figure |(7|) for two different mass splittings. The first plot is for small mass differences 
(AMA = 5 GeV and AMH C = 10 GeV), the second one for (relatively) larger splittings (AMA = 10 
GeV and AMH C = 50 GeV). In the second case, the cross section is too large to obtain an abundance 
consistent with WMAP, (cxv)wmap ~ pb. 

In the limit of small mass splittings and vanishing Xl = we have the right amount of dark matter 
provided Mh q ^,600 GeV. This regime corresponds to the narrow region around the diagonal in Figure 
(J6j) - The abundance increases for increasing Mh , but this can be somewhat compensated by playing 
with the mass splittings, which, as discussed above, tend to increase the cross section. For instance, 
for AMA = 5 GeV and AMH C = 10 GeV, we get the right relic abundance for M Ho > 800 GeV. 

There is however a limit to this trend. Indeed, the total annihilation cross section of a scalar 
particle, like our Hq, is constrained by unitarity to be smaller than 



a result which, in the context of dark matter relics from freeze-out, has been first put forward by Griest 
and Kamionkowski (32] ■ In the present model, increasing the mass splitting drives the annihilation cross 
section toward the strong quartic coupling regime. A similar result holds, for instance, for a weakly 
interacting massive neutrino candidate of dark matter. In this case, both the neutrino and its charged 
partner should be kept lighter than, say 1 TeV, since their mass comes from Yukawa couplings to the 
Higgs field. (See [33] and the discussion in [32]). In our case, the bulk of the mass of the components 
of the H2 doublet comes from the mass scale (X2, which is a priori arbitrary. However, the unitarity 
limit and WMAP data constraints, translate into the upper bound Mh < 130 TeV [32] . 

Unfortunately neither direct, nor indirect detection experiments are sensitive to the large mass 
region discussed in this section. Forthcoming direct detection experiments might do better, but as the 
plots show rather clearly, other forms of dark matter would then have to be introduced in order to 
explain the amount of dark matter that is currently observed. However, if the profile of dark matter 
in the galaxy is as assumed in the present paper (i.e. NFW), future gamma detectors (GLAST) will 
probe most of the parameter space considered in this section, keeping in mind that we had to finely 
tune the mass splittings to obtain the right abundance. 

4 Conclusions 

The dark matter candidate of the Inert Doublet Model stands fiercely by the neutralino. The lightest 
stable scalar is a weakly interacting massive particle with a rich, yet simple, phenomenology and it 
has a true potential for being constrained by existing and forthcoming experiments looking for dark 
matter. For the sake of comparison, we show in the right part of Figure J6]) a fair sample of models, 
both for the IDM and the MSSM. In the (MoM^DMh 2 ) plane, we clearly see the two regimes (low 
mass and high mass) of the IDM that may give rise to a relevant relic density (i.e near WMAP). The 
MSSM models have a more continuous behavior, with O(100 GeV) dark matter masses. For indirect 
detection, the IDM dark matter candidates have typically higher detection rates than the neutralino in 
SUSY models, especially at high mass. It is in particular interesting that the IDM can give the right 
relic abundance in a range of parameters which will be probed by GLAST. Direct detection is also 
promising but only at low masses, where the models are within reach of future ton-size experiments. 
The high mass models, however, have too small cross sections in the regions consistent with WMAP. 

The phenomenology of the IDM candidate for dark matter is intertwined with that of the Brout- 
Englert-Higgs particle, h. It could be of some interest to investigate the prospect for detection of the h 
and H2 components at the LHC, in the light of the constraints for dark matter discussed in the present 
paper. 
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Figure 5: From top to bottom: Relic density, gamma indirect detection and direct detection contours in 
the (M Ho ,ii 2 ) plane. Left: m h = 120 GeV. Right: m h = 200 GeV. 
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Figure 6: Left: same as Figure © for M Ho ,fi 2 G [700,1400] GeV. Right: Comparison with the MSSM. 
From top to bottom: Relic density, gamma indirect detection, direct detection, all as a function of the 
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Figure 7: (av) function of \i2 for Mh = 1000 GeV and small mass splittings AMAq = 5 GeV, 

AMH+ = 10 GeV for the plot on the left and larger mass splittings AMA = 10 GeV, AMH+ = 50 GeV 
for the plot on the right. Dashed line corresponds to annihilation into two Higgs, continuous lines to 
annihilation into two gauge bosons. 
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